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Glossary 

Ageing. A dynamical process involving a long-tailed waiting time distribution 
ijj{t) ~ ^a^-i-a (^0 < a < 1) does not possess a characteristic time scale t^{t)dt 
separating microscopic and macroscopic times. Instead, such a process exhibits distinct 
memory, such that the now-state of the system is strongly influenced by its state in the past. 
When the system is prepared at some time to < and the measurement starts at t = 0, this 
means that the result of the measurement depends on the time to even until rather long times. 



Anomalous diffusion. Under anomalous diffusion we understand here deviations of the linear 
time dependence {x^{t)) = 2Kt, of the mean squared displacement in absence of an external 
bias, in the form of a power-law: = 2K(t^. Here, K,^ is the anomalous diffusion 

constant of dimension cm^/sec''. In the range < C < 1, we deal with subdiffusion, whereas 
C > 1 describes super diffusion. Levy flights have a diverging mean squared displacement. 



Continuous time random walk. Continuous time random walk theory describes a ran- 
dom motion by assigning each jump a jump length x and a waiting time t elapsing in 
between two successive jumps, drawn from the two probability densities A(x) and i^it), 
respectively. The two densities A and ip fully specify the probability density function 
P(x, t) describing the random process. In Fourier-Laplace space, the propagator follows as 
Pik, u) = - i^{u)X{k)]. 



Fourier and Laplace transforms. Linear partial differential equations are often conveniently 
solved using integral transforms of the Laplace and Fourier type. Moreover, the definitions 
of the fractional operators used in the text correspond to Laplace and Fourier convolutions, 
such that these transformations also become useful there. The Laplace transform of a 
function f{t) is defined as 



f{u)^^{f{t);u}= I f{t)e-^'dt, 
while the Fourier transform of g{x) reads 

gik)^^{g{xy,k}= I g{x)e'^^dx. 

oo 



oo 



Note that we denote the transform of a function by explicit dependence on the respective 
variable. For m = and k = Q the Laplace and Fourier transform is but the average of the 
function f ot g, respectively. Tauberian theorems ascertain relations between the original 
function and its transform. For instance, the small u behaviour f{u) ~ 1 — [ut)° implies 
the long time scaling f{t) ~ r"/t^+". See Refs. (joj; 13) for details. 
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Fractional differintegration. The multiple derivative of an integer power is d'^x'^/dx"^ = 
n\/{n — m)!x"~™' for m > n. The result is zero if m > n. Replacing the factorials by 
the F-function, one can generalise this relation to qD^x^ = r(l + p)t^~'^ /r{l + p — q)- In 
particular, this includes the fractional differentiation of a constant, qD^I = x"'^/r(l — q), 
that does no longer vanish. Fractional differentiation was first mentioned by Leibniz in 
a letter to de I'Hospital in 1695. The Riemann-Liouville fractional operator used in the 
following is a similarly straightforward generalisation of the Cauchy multiple integral, 
followed by regular differentiation. 



Levy flight. A Levy flight is a special type of continuous time random walk. Its waiting 
time distribution is narrow for instance, Poissonian with ilj{t) = exp(— t/r), and the 
resulting dynamics therefore Markovian. The jump length distribution of a Levy flight is 
long-tailed: A(x) ~ cr'^|x|~^~'^, with < n < 2, such that no second moment exists. The 
resulting PDF is a Levy stable law with Fourier transform P{k,t) = exp (^—K^^H\k\^Y 



Levy stable laws. The generalised central limit theorem states that the properly normalised 
sum of independent, identically distributed random variables with finite variance converges 
to a Gaussian limit distribution. A generalisation of this theory exists for the case with 
infinite variance, namely, the generalised central limit theorem. The related distributions 
are the Levy stable laws, whose density in the simplest case have the characteristic function 
p{k) = J^^p{x) exp{ikx)dx (Fourier transform) of the form p{k) = exp (— c|/c|'^), with 
< /i < 2. In direct space, this corresponds to a power-law asymptotic p{x) ^ c|x|~^~'^. In 
the limit /i = 2, the universal Gaussian distribution is recovered. 



Levy walks. In contrast to Levy flights. Levy walks possess a finite mean squared displace- 
ment, albeit having a broad jump length distribution. This is possible by the introduction 
of a time penalty for long jumps through a coupling X{x)p{t\x) between waiting times 
and jump lengths, such that long jumps involve a longer waiting time. For instance, a 
(5-coupling of the form ^X{x)6{\x\—vt) is often chosen, such that v plays the role of a velocity. 



Strange kinetics. Often, deviations from exponential relaxation patterns exp(— t/r) and 
regular Brownian diffusion are observed. Instead, non-exponential relaxation, for instance, 
of the stretched exponential form exp(— [t/r]") (0 < a < 1) or of the inverse power-law 
form t~^, are observed, or anomalous diffusion behaviour is found. 



Weak ergodicity breaking. In a system with a broadly distributed waiting time with diverging 
characteristic time scale, a particle can get stuck at a certain position for a long time. For 
instance, for a waiting time distribution of the form tlj{t) ~ r°t~^~", the probability of not 
moving until time t scales as T°'t~°', i.e., decays very slowly. The probability for not moving 
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is therefore appreciable even for long times. A particle governed by such a ip{t) even at 
stationarity does not equally explore a different domains of phase space. 



I. DEFINITION AND IMPORTANCE OF ANOMALOUS DIFFUSION 



Classical Brownian motion characterised by a mean squared displacement 

{x\t)) = 2Kt, 



(1) 



growing linear in time in absence of an external bias is the paradigm for random motion. 
Note that we restrict our discussion to one dimension. It quantifies the jittery motion of 
coal dust particles observed by Dutchman Jan Ingenhousz in 1785 (1l4l ). the zigzagging of 
pollen grain in solution reported by Robert Brown in 1827 (jZlh . and possibly the dance 
of dust particles in the beam of sunlight in a stairwell so beautifully embalmed in the 
famed poem by Lucretius (0). In general, Brownian motion occurs in simple, sufficiently 
homogeneous systems such as simple liquids or gasses. In the continuum limit, Brownian 
motion is governed by the diffusion equation 



dt 



(2) 



for the probability density function (PDF) P{x,t) describing the probability P{x,t)dx to 
find the particle at a position in the interval x, . . . ,x + dx aX time t. For a point-like initial 
condition P{x, 0) = S{x), the solution becomes the celebrated Gaussian PDF 



P{x,t) 



exp 



X 



(3) 



The diffusion constant K fulfils the Einstein-Stokes relation K = ksT / {mri), where ksT is 
the Boltzmann energy at temperature T, m is the mass of the test particle, and rj the friction 
coefficient. This relation between microscopic and macroscopic quantities was used for the 
determination of the Avogadro number by Perrin (1971 ). Examples of Perrin's recorded ran- 
dom walks and the jump length distribution constructed from the data for a time increment 
of 30 sec are shown in Fig. [TJ Fig. [2] displays time traces of the Brownian motion of a small 
mirror in air that lead to an amazingly accurate determination of the Avogadro number by 
Keppler, for two different ambient pressures (1651 ). 

However, i n ma ny systems deviations from the linear behaviour ([1]) are observed (0; 0; 



17l : Il8l : l29l : l68l : llOQl ). These deviations can assume the power- law form 



{x\t)) 



(4) 



For < C < 1; we observe subdiffusion. Prominent examples of subdiffusion include charge 
carrier tr ansp ort in amorphous semiconductors (9^; 1031) . diffusion of chemicals in subs urfac e 
aquifers ( 104 ). the motion of beads in actin gels (IiIQI ). motion in chaotic maps (j3 ^ 120l ). 
or the subdiffusion of biomacromolecules in cells (j59l): just to name a few [compare (jl8l) for 
more details]. Subdiffusion of this type is characterised by a long-tailed waiting time PDF 



-l-Q 



corresponding to the time-fractional diffusion equation 



dP{x,t) 

Ft 



KaoD 



1-a 



dx"^ 



P{x,t), 



(5) 



Fig. 6. 
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FIG. 1 Random walk trajectories of putty particles in water recorded by Perrin (jOSl ). Left: three 
designs obtained by tracing a small grain of putty at intervals of 30sec. Right: the starting point 
of each motion event is shifted to the origin. The figure illustrates the continuum approach of the 
jump length distribution if only a large number of jumps is considered. 



with the fractional Riemann-Liouville operator (1201 : l2ll : |2J) 



Dl-^P{x,t) 



d 



T{a) dt Jo {t - t 



(6) 



From the latter definition, it becomes apparent that subdiffusion corresponds to a slowly 
decaying memory integral in the dynamical equation for P{x,t). 

The lack of a characteristic time scale (i.e., the divergence of t%l){t)dt) of the waiting 
time PDF %l){t) ^ no longer permits to distinguish microscopic and macroscopic events. 
This causes that the diffusing particle can get stuck at a certain position for very long 
times, quantified by the sticking probability (pit) = J^°° ip{t')dt' of not moving. The PDF 
P(x, t) exhibits characteristic cusps at the location where the particle was initially released 
and, in the presence of an external drift, growing asymmetry. Moreover, a process whose 
measurement be gins at t = depends on the preparation time at some prior time, the 
so-called ageing (1361 : l95l ). Another effect in this context is that a particle no longer evenly 
distributes in a certain space, such that despite the existence of stationary states a weak 
ergodicity breaking occurs (38; [H; 77). 

Contrasting a subdiffusing particle are Levy flights (LFs), that are based on a waiting 
time PDF with finite characteristic time but a jump length distribution \{x) ~ 
(0 < /i < 2) with diverging variance \x\'^\{x)dx. The diffusion equation for an LF 
becomes generalised to the space-fractional diffusion equation 



(7) 
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FIG. 2 Erroneous behaviour of Brownian motion made visible in a high-precision measurement. 
Data from an Edelmann recorder, obtained by Kappler in 1931 (|65l ). Kappler monitored the Brow- 
nian motion of a small mirror (surface approx. Imm^), suspended from a fine quartz thread (several 
cm long and some /xm thick). The mean squared of the torsional displacement, (p"^, follows the 
relation Dip"^ = ksT, where D is the directional force of the suspension (|65l ). The facsimiles show 
three different realisations. From his data, Kappler obtained the Avogadro-Loschmidt number 
Nl = 60.59 • 10^2 ± 1%, to a remarkable accuracy. 



where the fractional Riesz-Weyl operator is defined through its Fourier transform 
^ {d'^ / d\x\'^ P {x , t)} = — |A;|'*P(/c, t). In Fourier space, one therefore recovers immediately 
the characteristic function P{k,t) = exp (^—K^^^\k\^t^ of a symmetric Levy stable law. As 
can be seen from Eq. LFs are Markovian processes. However, their trajectory now has a 
fractal dimension df = ^ and is characterised by local search changing with long excursions, 
see Fig. |3l This property has be en sh own to be a better search strate gy th an Brownian mo- 
tion, as it oversamples less (76; 115 ). I n fac t, animals like albatross ( 114l ). spider monkeys 
fflOoh . jackals (HH), and even plankton (11131 ) and bacteria (1751 ) were claimed to follow Levy 
search strategies. LFs also occur in diffusion in energy space (Il23l ). or in optical lattices 
(1661 ). Due to the clustering nature of their trajectory, LFs also exhibit a form of ergodicity 
breaking (jTol ). 



II. INTRODUCTION 



For sums of independent, identically distributed random variables with proper normali- 
sation to the sample size, the generalised central limit theorem guarantees the convergence 
of the associated probability density to a Levy stable density even though the variance of 
these random variables diverges (j^; 11; 13; [la; 26). Well-known examples for Levy stable 



FIG. 3 Trajectories of Brownian motion (left) and a Levy flight of index /_f = 1.5 (right), both 
with the same number (~ 7000) of steps. The long sojom'ns and clustering appearance of the Levy 
flight are distinct. 



densities are the one-sided (defined for x > 0) Levy-Smirnov distribution 



/iA-i/2(^) = \/^e-p(-^)' (8) 



related to the first passage time density of a Gaussian random walk process of passing the 
origin (see below), and the Cauchy (or Lorentz) distribution 

fM = -tA-^v (9) 

vr (1 + x^) 

In general, an Levy stable density is defined through its characteristic function of the PDF 

^(z) ^ ^ {fix)} = / f,A^)e^'''dx (10) 



where 

\og(p{z) = — |2;|^exp |z-;^sign(2;) /> , (11) 



2 

for yU 7^ 1. Here, the skewness (or asymmetry) parameter /5 is restricted to the following 
region: 

l/?l<l?' (12) 

- I 2 - /i, if 1 < ;U < 2. ^ ^ 

For (3 = 0, the corresponding Levy stable density is symmetric around x = 0, while for 
j3 = — /i and < /i < 1, it is one-sided. In general, an Levy stable density follows the 
power-law asymptotic behaviour 



A 



ff^A^) - Tzf&.^ < 2, (13) 



\x 



i+fj. 



with Afj^^/s being a constant, such that for all Levy stable densities with fi < 2 the variance 
diverges 

(x^) = oo. (14) 

Conversely, all fractional moments (|a;|'^) < oo for all < 5 < < 2. From above definitions 
it is obvious that the Levy stable density /2,o corresponds to the Gaussian normal distribution 



f2,o{x) = \l^ exp ( --X ] (15) 
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possessing finite moments of any order. In tliis limit, tlie generalised central limit theorem 
coincides with the more traditional, and universal, central limit theorem. 

Brownian motion has traditionally been employed as the dominant model of choice for 
random noise in continuous-time systems, due to its remarkable statistical properties and 
its amenability to mathematical analysis. However, Brownian motion is just a single ex- 
ample of the Levy family. Moreover, it is a very special and somewhat misrepresenting 
member of this family. Amongst the Levy family, the Brownian member is the only motion 
with continuous sample-paths. All other motions have discontinuous trajectories, exhibiting 
jumps. Moreover, the Levy family is characterized by selfsimilar motions. Brownian motion 
is the only selfsimilar Levy motion possessing finite variance, while all other selfsimilar Levy 
motions have an infinite variance. 

Random processes whose spatial coordinate x or clock time t are distributed according 
to an Levy stable density exhibit anomalies, that is, no longer follow the laws of Brownian 
motion. Consider a continuous time random walk pr oces s defined in terms of the jump length 



and waiting time distributions X{x) and tlj{t) (1961 : Il03l ). Each jump event of this random 
walk, that is, is characterised by a jump length x drawn from the distribution A, and the 
time t between two jump events is distributed according to ip. (Note that an individual 
jump is supposed to occur spontaneously.) In absence of an external bias, continuous time 
random walk theory connects X{x) and ip{t) with the probability distribution P{x,t)dx to 
find the random walker at a position in the interval (x, x + dx) at time t. In Fourier-Laplace 



space, P{k, u) = ^ {=Sf {P(x, t); t — >• m} ; x — >• /c}, this relation reads (1671 ) 



Pft^)- ^ (16) 

u 1 — X{k)ip[u) 

where J^{f{t)} = exp(—ut)f{t)dt. We here neglect potential complications due to 
ageing effects. The following cases can be distinguished: 

(i) A(x) is Gaussian with variance cr^ and ip{t) = 6{t — t). Then, to leading order in k"^ 
and u, respectively, one obtains X{k) ^ 1 — a^k"^ and ip{u) ~ 1 — ut. From relation (fT6|) 
one recovers the Gaussian probability density P{x,t) = ^/T/l4^^Kt) exp{—x'^ / (AKt)} with 
diffusion constant K = a'^/r. The corresponding mean squared displacement grows linearly 
with time, see Eq. ([T]). This case corresponds to the continuum limit of regular Brownian 
motion. Note that here and in the following, we restrict the discussion to one dimension. 

(ii) Assume X{x) still to be Gaussian, while for the waiting time distribution ip{t) we 
choose a one-sided Levy stable density with stable index < a < 1. Consequently, 'ipi^u) ~ 
1 — (ut)", and the characteristic waiting time tip{t)dt diverges. Due to this lack of a time 
scale separating microscopic (single jump events) and macroscopic (on the level of P jx, t) ) 



scales, P{x,t) is no more Gaussian, but given by a more complex if-function (llTl : Il8l : llO 



In Fourier space, however, one finds the quite simple analytical form (llTI ) 



P{k, t) = K i-K^kH'^) = Y: rn V \ (17) 

^ ^ ^ 1(1 + Q^^j 

in terms of the Mittag-Leffler function. This generalised relaxation function of the Fourier 
modes turns from an initial stretched exponential (KWW) behaviour 

P{k,t) ~ 1- K^Pr/Til + a) ^ exp {-Kc^Pr /T{1 + a)} (18) 
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to a terminal power-law behaviour (jlTl ) 



P{k,t) ~ (^K„rrr(l - a)j . (19) 

In the limit a — > 1, it reduces to the traditional exponential P{k,t) = exp{—Kk'^t) with 
finite characteristic waiting time. Also the mean squared displacement changes from its 
linear to the power-law time dependence 

{x\t)) = (20) 

with Kct = c'^/tq. This is the case of subdiff usion . We note that in x, t space the dynamical 



equation is the fractional diffusion equation (11051 ). In the presence of an external potential 



it generalises to the time- fractional Fokker-Planck equation (llTl : 1 181 : l83l ). see also below. 

(iii) Finally, take ipit) = S{t — T) sharply peaked, but X{x) of Levy stable form with index 
< fi < 2. The resulting process is Markovian, but with diverging variance. It can be 
shown that the fractional moments scale like fjo^ ) 

{\x{t)\') cx {K'^^kf , (21) 

were 

j^M = ^M/^. The upper index is chosen to distinguish from the subdiffusion con- 
stant K^. Note that the dimension of K^^^ is cm'^/sec. From Eq. ( |T6i) one can immediately 
obtain the Fourier image of the associated probability density function, 

P{k,t) =expi^- K'^^'^lkl^y (22) 

From Eq. fllip this is but a symmetric Levy stable density with stable index /i, and this 
type of random walk process is most aptly coined a Levy flight. A Levy flight manifestly 
has regular exponential mode relaxation and is in fact Markovian. However, the modes 
in position space are no more sharply localised like in the Gaussian or subdiffusive case. 
Instead, individual modes bear the hallmark of an Levy stable density, that is, the diverging 
variance. We will see below how the presence of steeper than harmonic external potentials 
causes a flnite variance of the Levy flight, although a power- law form of the probability 
density remains. 

In the remainder of this paper, we deal with the physical and mathematical properties 
of Levy flights and subdiffusion. While mostly we will be concerned with the overdamped 
case, in the last section we will address the dynamics in velocity space in the presence of 
Levy noise, in particular, the question of the diverging variance of Levy flights. 



III. LEVY FLIGHTS 

A. Underlying random walk process 

To derive the dynamic equation of a Levy flight in the presence of an external force fleld 
F{x) = —dV{x)/dx, we pursue two different routes. One starts with a generalised version of 
the continuous time random walk, compare Ref. (1841 ) for details; a slightly different derivation 
is presented in Ref. fl37l ). 

To include the local asymmetry of the jump length distribution due to the force 
field F{x), we introduce (1841 : \91\ ] the generalised transfer kernel A{x,x') = \{x — 
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x') [A{x')Q{x — x') + B{x')Q{x' — x)] (and therefore A(x, x') = A{x'; x — x')). As in stan- 
dard random walk theory (compare (|28|)). the coefficients A{x) and B{x) define the local 
asymmetry for jumping left and right, depending on the value of F{x). Here, G(x) is the 
Heaviside jump function. With the normalisation J A{x', A)dA = 1, the fractional Fokker- 
Planck equation (FFPE) ensues (0): 



d_ 

dt 



P{x,t) 




dx mrj 



(23) 



Remarkably, the presence of the Levy stable A (a;) only affects the diffusion term, while the 



drift term remains unchanged (l56l : l57l : l84j ). The fractional spatial derivative represents an 



integrodifferential operator defined through 



9^ 
^\x\^ 



■P{x,t) 



2cos(7r/i/2)r(2-/i)(9x2 



P{x\t) 



\x 



X' 



(24) 



for 1 < /i < 2, and a similar form for < /^i < 1 (1211 : |2J; |45| ). In Fourier space, for all 
< /X < 2 the simple relation 



d\xY 



■P{x,t) 



\k\^P{k,t) 



(25) 



holds. In the Gaussian limit /i = 2, all relations above reduce to the familiar second-order 
derivatives in x and thus the corresponding P{x, t) is governed by the standard Fokker- 
Planck equation. 
The FFPE 



can also be derived from the Langevin equation (j45|; l56|; l57|; |64| ) 



dx(t) 
dt 



1 dV{x) 



dx 



(26) 



driven by white Levy stable noise ^f,{t), defined through L{At) = jl^"^' i^{t') dt' be- 
ing a symmetric Levy stable density of index /x with characteristic function p(/c. At) = 
exp (— K ^\k\^ At) for < /i < 2. As with standard Langevin equations, K^^'^ denotes the 
noise strength, m is the mass of the diffusing (test) particle, and 7 is the friction constant 
characterising the dissipative interaction with the bath of surrounding particles. 

A subtl e poi nt about the FFPE (!23|) is that it does not uniquely define the underlying 
trajectory (jlOSi ): however, starting from our definition of the process in terms of the stable 
jump length distribution A(x) ~ or its generalised pendant A{x,x'), the FFPE ([23]) 

truly represents a Levy fiight in the presence of the force F{x). This poses certain difficulties 
when non-trivial boundary conditions are involved, as shown below. 



B. Propagator and symmetries 



In absence of an external force, F{x) = 0, the exact solution of the FFPE is readily 
obtained as the Levy stable density P(k,t) = exp (— i^^^^ |A;|'^t) in Fourier space. Back- 
transformed to position space, an analytical solution is given in terms of the Fox if-function 



P{x,t) 



^\x\ 



1,1 
2,2 



\x\ 



{K(f'H) 



;i,i//i), (1,1/2) 

;i,i), (1,1/2) 



(27) 
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Cauchy, t=0.1 
Cauchy, t=0.5 
Gauss, t=0.1 
Gauss, t=0.5 



0.1 



0.01 




FIG. 4 Cauchy distribution (/x = 1) for two times in comparison to the Gaussian (/i = 2). We chose 
K^^") = K = 1. Note that the Cauchy distribution is narrower at the origin, and after crossing the 
Gaussian falls off in the much slower power-law fashion. 



from which the series expansion 



P{x,t) = 
derives. For fi 



E 



r(|i + 



\X\ 



/x(KMt)i^^F([l + z/]/2)F(l-[l + z/]/2) u\ \{K(^'Hy/^' 
= 1, the propagator reduces to the Cauchy Levy stable density 

P{x,t) 



(28) 



1 



(29) 



We plot the time evolution of P{x,t) for the Cauchy case /i = 1 in Fig. |l]in comparison to 
the limiting Gaussian case /i = 2. 

Due to the point symmetry of the FFPE (l23l) for F{x) = 0, the propagator P(x, t) is 
invariant under change of sign, and it is monomodal, i.e., it has its global maximum at x = 0, 
the point where the initial distribution P(x, 0) = S{x) was launched at t = 0. The latter 
property is lost in the case of strongly confined Levy flights discussed below. Due to their 
Markovian character. Levy flights also possess a Galilei invariance (17; sl). Thus, under the 
influence of a constant force field F{x) = Fq, the solution of the FFPE can be expressed in 
terms of the force-free solution by introducing the wave variable x — Fgi, to obtain 



Ppo(x,t) 



Pq { X , t 

7717 



(30) 



This result follows from the FFPE fl23p . that in Fourier domain becomes f l64l ) 



(31) 
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0.1 



0.01 




Cauchy, t=0.5 
Cauchy, t=2.0 
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FIG. 5 Cauchy distribution with K^-^^ = 1 advected along a field FQ/{mj) = 2, for different times. 



with solution 



P{k, t) = exp 



m7 



(32) 



By the translation theorem of the Fourier transform, Eq. fl30l) yields. We show an example 
of the drift superimposed to the dispersional spreading of the propagator in Fig. [51 



C. Presence of external potentials 

1. Harmonic potential 

In an harmonic potential V{x) = ^Xx^, an exact form for the characteristic function can 
be found. Thus, from the corresponding FFPE in Fourier space. 



m7 dk 

by the method of characteristics one obtains 



P{k, t) = exp 



[1 



(33) 



(34) 



for an initially central (5-peak, P(a;, 0) = 5{x) (1641 ). This is but the characteristic function 
of an Levy stable density with time- varying width. For short times, 1 — exp(— /iAt/[m7]) ~ 
jj.Xt/lm'y] grows linearly in time, such that P{k,t) ~ exp (^—K^^^\k\^tj as for a free Levy 
flight. At long times, the stationary solution defined through 
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is reached. Interestingly, it has the same stable index as the driving Levy noise. By 
separation of variables, a summation formula for P{x, t) can be obtained similarly to the 
solution of the Ornstein-Uhlenbeck process in the presence of white Gaussian noise, however, 
with the Hermite polynomials replaced by iZ-functions (jG^ ). 

We note that in the Gaussian limit /i = 2, the stationary solution by necessity has to 
match the Boltzmann distribution corresponding to exp (— /crT/c^/PA]). This requires that 
the Einstein-Stokes relation K = ksT /{m'^) is fulfilled (1271 ). One might therefore speculate 
whether for a system driven by external Levy noise a generalised Einstein- Stokes relation 
should hold, as was established for the subdiffusive case (jlTI : l83l ). As will be shown now, 
in steeper than harmonic external potentials, the stationary form of t) even leaves the 
basin of attraction of Levy stable densities. 



2. Steeper than harmonic potentials 

To investigate the behaviour of Levy flights in potentials, that are steeper than the 
harmonic case considered above, we introduce the non-linear oscillator potential 

V{x) = + ^-x\ (36) 

that can be viewed as a next order approximation to a general confining, symmetric potential. 
It turns out that the resulting process differs from above findings if a suitable choice of the 
ratio a/h is made. For simplicity, we introduce dimensionless variables through 



where 



arriving at the FFPE 



x^x/xq] t^t/to] a ^ ato / (mj) , (37) 



|P(,,<) + |,|.= (,_^_,,|.)P(,,<). (39) 



Consider first the simplest case of a quartic oscillator with a = in the presence of Cauchy 
noise {fi = 1). In this limit, the stationary solution can be obtained exactly, yielding the 
expression 

1 

n 1 — x'^ + x"^ 

plotted in Fig. [6l Two distinct new features in comparison to the free Levy flight, and the 
Levy flight in an harmonic potential: (1) Instead of the maximum at x = 0, one observes 
two maxima positioned at 

Xm = ±Vlj2; (41) 
at X = 0, we find a local minimum. (2) There occurs a power-law asymptote 

W ~ ^ (42) 
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for X ^ 1; consequently, this stationary solution no longer represents an Levy stable density, 
and the associated mean squared displacement is finite, (x^) < oo. 

A more detailed analysis of Eq. fl59|l reveals (jH; IIHI ). that (i) the bimodality of P{x,t) 
occurs only if the amplitude of the harmonic term, a, is below a critical value a^, (ii) for 
general yU, the asymptotic behaviour is Pgt and there 

exists a finite bifurcation time tc at which the initially monomodal form of t) acquires 
a zero curvature at x = 0, before settling in the terminal bimodal form. 

Interestingly, in the more general power-law behaviour 

V{x) = ^, (43) 

the turnover from monomodal to bimodal form of P{x,t) occurs exactly when c > 2. The 
harmonic potential is therefore a limiting case when the solution of the FFPE still belongs 
to the class of Levy stable densities and follows the generalised central limit theorem. This is 
broken in a superharmonic (steeper than harmonic) potential. The corresponding bifurcation 



time tc is finite for all c > 2 (l45l ) . An additional effect appears when c > 4: there exists a 
transient trimodal state when the relaxing S{x)-peak overlaps with the forming humps at 
X = ±Xm,. At the same time, the variance is finite, if only c > 4 — fi, following from the 
asymptotic stationary solution 

W ~ (44) 



TTlX 



Details of the asymptotic behaviour and the bifurcations can be found in Refs. (1431 : l45l ). 



From a reverse engineering point of view. Levy fiights in confining potentials are studied in 
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D. First passage and first arrival of Levy flights 

One might naively expect that a jump process of Levy type, whose variance diverges 
(unless confined in a steep potential) may lead to ambiguities when boundary conditions 
are introduced, such as an absorbing boundary at finite x. Indeed, it is conceivable that for 
a jump process with extremely long jumps, it becomes ambiguous how to properly define 
the boundary condition: should the test particle be absorbed when it arrives exactly at the 
boundary, or when it crosses it anyplace during a non-local jump? 

This question is trivial in the case of a narrow jump length distribution: all steps are small, 
and the particle cannot jump across a point (in the continuum limit considered herein). For 
such processes, one enforces a Cauchy boundary condition P(0, t) = at the point x = 
of the absorbing boundary, removing the particle once it hits the barrier after starting at 
Xq, where the dynamics is governed by Eq. (125]) with F{x) = 0. Its solution can easily 
be obtained by standard methods, for instance, the method of images. This is completely 
equivalent to considering the first arrival to the point x = 0, expressed in terms of the 
diffusion equation with sink term: 

-{x, t) = K—^{x, t) - Pfa(t)5(x), (45) 



dt ' dx'^ 

defined such that P(0, t) = 0. Note that the quantity ^ is no longer a probability density, 
as probability decays to zero; for this reason, we use the notation J^. From Eq. (H5l) by 
integration we obtain the survival probability 

(46) 

with ^(0) = 1 and limf^oo -^(^) =0. Then, the first arrival density becomes 

pUt) = (47) 

Eq. (H5!) can be solved by standard methods (determining the homogeneous and inhomoge- 
neous solutions). It is then possible to express ^{x,t) in terms of the propagator P{x,t), 
the solution of Eq. (1231) with F{x) = with the same initial condition, P{x, 0) = 6{x — Xq) 
and natural boundary conditions. One obtains 

P(0,t)= / pUr)P{xo,t-T)dT, (48) 
Jo 

such that the first arrival density corresponds to the waiting time distribution to jump from 
Xq to (or, vice versa, since the problem is symmetric). In Laplace space, this relation takes 
on the simple algebraic form pf^{u) = P(0, u)/P{xo, u). Both methods the explicit boundary 
value problem and the first arrival problem for Gaussian processes produce the well-known 
first passage (or arrival) density of Levy-Smirnov type ([8]), 

pit) = pi,{t) = ^ exp (-^] ~ , ^° , (49) 

^ ^ ^ ^ 74^ V ^KtJ ^ ^ 

with the asymptotic power- law decay p{t) ~ t~'^^'^, such that no mean first passage time 



exists (jl3l : l23l ) . 
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FIG. 7 First passage density for various stable indices /i and initial position xq = 10.0 away from 
the absorbing boundary. Again, the universal ~ t~^^'^ scaling is distinct. Note that we plot tp{t) 
in the y-axis. 



Long-tailed jump length distributions of Levy stable form, however, endow the test parti- 
cle with the possibility to jump across a certain point repeatedly. The first arrival necessarily 
becomes less efficient. Indeed, as shown in Ref. (1441 ). the Gaussian result (H^ is generalised 
to 



as t 



oo 



(50) 



with C(/i) = /ir(2 - /i)r(2 - l//i) sin(7r[2 - /i]/2) sin^in/ - 1]), and 1 < /i > 2 (j4J). 
The long-time decay ~ ^-2+1/^* -g giQ-^g^ than in (jUj). 

One might naively assume that the first passage problem (the particle is removed once 
it crosses the boundary) for Levy flights should be more efficient, that is, the first passage 
density p{t) should decay quicker, than for a narrow jump length distribution. However, as 
we have a symmetric jump length distribution A (a;), the long outliers characteristic for these 
Levy flights can occur both toward and away from the absorbing barrier. From this point of 
view it is not totally surprising to see the simulations result in Fig. [TJ that clearly indicate 
a universal asymptotic decay ~ exactly as for the Gaussian case. 

In fact, for all Markovian pr oces s es w ith a symmetric jump length distribution, the 
Sparre Andersen theorem ([^; [12; 110 : llll ) proves without knowing any details about X{x) 
the asymptotic behaviour of the first passage time density universally follows p{t) ~ t""^/^. 
The details of the specific form of A(x) only enter the prefactor, and the pre- asymptotic 
behaviour. A special case of the Sparre Andersen theorem was proved in Ref. (jTol ) when the 
particle is released at Xq = at time t = 0, and after the first jump an absor bing boundary 
is installed at x = 0. This latter case was simulated extensively in Ref. (jl2ll ). From a 
fractional diffusion equation point of view, it was shown in Ref. (jH) that the fractional 
operator d'^/d\x\'^ needs to be modified, to account for the fact that ^(x, t) = beyond the 
absorbing boundary, such that long-range correlations are present exclusively for all x in the 
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semi-axis containing xq. The fractional diffusion equation in the presence of the absorbing 
boundary therefore has to be modified to (1441 ) 



d K 92 



\x — X 



dx'. 



(51) 



where k = 2r(2 — /i) |cos(7r/i/2)|, such that the first term on the right hand side no longer 
represents a Fourier convolution. An approximate solution with Cauchy boundary condition 
reveals p(u) ~ 1 — cm^/^, where c is a constant, indeed leading to the Sparre Andersen 
behaviour p(t) ~ 

This also demonstrates that the method of images no longer applies when Levy flights 
are considered, for the images solution 



i(a;, t) = P{x - xo, t) - P{x + xq, t) 



(52) 



would be governed by the full fractional diffusion equation, and not Eq. (I5T]) . and the result 
for the first passage density, p{t) ~ would decay faster than the Sparre Andersen 

universal behaviour. A detailed discussion of the applicability of the method of images is 
given in terms of a subordination argument in Ref. (jlOSi ). We emphasise that this subtle 



failure of the method of images has been overlooked in literature previously (jl9l : l58l ). and 



care should therefore be taken when working with results based on such derivations. We 
also note that the method of images works in cases of subdiffusion, as the step length is 
narrow (l90l ). 



E. Leapover properties of Levy flights 

The statistics of first passage times is a classical concept to quantify processes in which it 
is of interest when the dynamic variable crosses a certain threshold value for the first time. 
For processes with broad jump length distributions, another quantity is of interest, namely, 
the statistics of the first passage leapovers, that is, the distance the random walker overshoots 
the threshold value x = d in a single jump (see Fig. [HD. Surprisingly, for symmetric LFs with 
jump length distribution A(x) ~ |x|~^~^ (0 < /U < 2) the distribution of leapover lengths 
across Xq is distributed like pi{t) ~ i.e., it is much broader than the original jump 

length distribution. In contrast, for one-sided LFs jumps the scaling of pi{€) bears the same 
index Information on the leapover behaviour is important to the understanding of how 
far search processes of animals for food or of proteins for their specific binding site along 
DNA overshoot their target, or to define better stock market strategies determining when 
to buy or sell a certain share instead of a given threshold price. 

Using a general theorem for the first passage times and leapovers for homogeneous pro- 
cesses with independent increments, leapover properties for completely symmetric and fully 



asymmetric (one-sided) LFs are derived in (1711 ). the general case is considered in (172|). The 
basic results are as follows. For the completely symmetric LF with index /x, the distribution 
of first passage leapover lengths I for a particle originally released a distance d away from 
the boundary reads (in scaled units) 



(53) 
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-X 



x = 

FIG. 8 Schematic of the leapover problem: the random walker starts at x = and after a number 
of jumps crosses the point x = d, overshooting it by a distance £. For narrow jump length 
distributions, each jump is so small that crossing of the point d is equal to arriving at this point. 



The validity of this result is confirmed by extensive simulations, for more details refer to 



( l70l : l71l : 172). Note that pi is normalised. In the limit fx ^ 2, pi tends to zero if £ 7^ and 
to infinity at ^ = corresponding to the absence of leapovers in the Gaussian continuum 
limit. However, for < yU < 2 the leapover PDF follows an asymptotic power-law with index 
/i/2, and is thus broader than the original jump length PDF A(x) with index fi. This is a 
remarkable finding: while A for 1 < /i < 2 has a finite characteristic length this always 

diverges for pi{i) irrespective of fi. 



In contrast, the result for completely asymmetric LFs has the form (1521 : l7ll ) 



sin(7r/i) 



TT 



(54) 



leading to the leapover PDF 



sin(7ryu) 



TT 



it^id + i) 



(55) 



which corresponds to the result obtained in Refs. (5^; 54) from a different method. Thus, for 
the one-sided LF, the scaling of the leapover is exactly the same as for the jump length dis- 
tribution, namely, with exponent /i. Again, this result compares favourably with simulations 



F. Kramers problem for Levy flights 

Many physical and chemical problems are related to the thermal fluctuations driven 
crossing of an energetic barrier, such as dissociation of molecules, nucleation processes, or 
the escape from an external, confining potential of finite height (jl2l ). A particular example 
of barrier crossing in a double-well potential driven by Levy noise was proposed for a long 
time series of paleoclimatic data (jsol ). Further cases where the crossing of a potential barrier 
driven by Levy noise is of interest is in the theory of plasma devices ( l42l ) , among others (jlSl ) . 

To investigate the detailed behaviour of barrier crossing under the influence of external 
Levy noise, we choose the rather generic double well shape 



Vix) = + \x\ 



(56) 
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t 

FIG. 9 Probability density function p{t) of barrier crossing times for fi = 1.0 and D = 10"^'^ 
0.00316. The dashed hne is a fit to equation 1^ with MOT = 1057.8 ± 17.7. 



Integrating the Langevin equation (l26l) with white Levy noise, we find an exponential decay 
of the survival density in the initial well: 

p(t) = ^exp(^-^^, (57) 

as demonstrated in Fig. [9l Levy flight processes being Markovian, this is not surprising, 



since the mode relaxation is exponential (jlTl : Il8l ). More interesting is the question how 



the mean escape time behaves as function of the characteristic noise parameters D and 
fi. While in the regular Kramers problem with Gaussian driving noise the Arrhenius-type 
activation Tc = Cexp{h/D) is followed, where h is the barrier height, and the prefactor C 
includes details of the potential, in the case of Levy noise, a power-law form 



Tc{^^^D) = —^^ (58) 



was assumed fl46l). Detailed investigations (1481 ) show that the scaling exponent uiu) = 1 



for all /i strictly smaller than 2. As already proposed in Ref. (15 ll ) and derived in (1631 ) in 
a somewhat different model, this means that, apart from a prefactor, the Levy flight is 
insensitive to the external potential for the barrier crossing, as confirmed by simulations 
(48l). Note that in comparison to Ref. (H), also values of /i in the range (0, 1) are included. 
For large values of deviations from the scaling are observed: eventually it will only take a 
single jump to cross the barrier when D oo. Detailed studies show indeed that eventually 
the unit time step is reached, i.e., Tc — 1. 
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G. More on the "pathology" 



Despite their mathematical foundation due to the generahsed central limit theorem and 
their broad use in the sciences and beyond as description for statistical quantities, and 
despite the existence of systems (for instance, the diffusion on a po lyme r in chemical space 



106|)), the divergence of 



mediated by jumps where the polymer loops back on itself fl40l : \7t 
the fluctuations of Levy processes is sometimes considered a pathology. This was already 
put forward by West and Seshadri (11161 ). who pointed out that a Levy flight in velocity 
space would be equivalent to a diverging kinetic energy. Here, we show that higher order 
dissipation effects lead to natural cutoffs in Levy processes. 

At higher velocities the friction experienced by a moving body starts to depend on the 
velocity itself (0). Such non-linear friction is known from the classical Riccati equation 
Mdv{t) / dt = Mq — Kv{tY for the fall of a particle of mass M in a gravitational field with 
acceleration o (|8|), or autonomous oscillatory systems with a friction that is non- linear in 
the velocity (|2|; y). The occurrence of a non-constant friction coefficient 7(^) leading to a 
non-linear dissipative force —'^{y)V was highlighted in Klimontovich's theory of non-linear 
Brownian motion (jlSl ). It is therefore natural that higher order, non-linear friction terms 
also occur in the case of Levy processes. 

We consider the velocity-dependent dissipative non-linear form (necessarily an even func- 
tion) (0) 



l{V) = 70 + + ... + 12nV 



2N 



l2N > 



(59) 



for the friction coefficient of the Levy flight in velocity space as governed by the Langevin 
equation 

dV{t) + -f{V)V{t)dt = dL{t) (60) 

with the constant friction 70 = 7(0). L{t) is the /i-stable Levy noise defined in terms 
of a characteristic function p*(a;,t) = ^{p{L,t)} = f^^p{L,t) exp (iuL) dL of the form 



p*{uj,t) = exp {—D\uj\H) (jl6l : l25l : l26l ). where D of dimension cm^/sec is the generalised 



diffusion constant. This is equivalent to the fractional Fokker Planck equation (1171 : Il8l : |4 



45|; l56|; l57|) 



dP{V, t) 



±{V,(V)P)+D 



(61) 



dt dV\ " ' J ^\V\^' 

As we showed in Sec. 1III.C.21 by the example of the Levy flights in position space, the 
presence of the first higher order correction, 72^^ in the friction coefficient 'y{V) rectifies the 
Levy motion such that the asymptotic power-law becomes steeper and the variance finite. 
When even higher order corrections are taken into consideration, also higher order moments 
become finite. We show an example in Fig. [10] for the second moment. 

The effect on the velocity distribution of the process defined by Eqs. (I6OI) and (|6T]) for 
higher order corrections are demonstrated in Fig. [11] for the stationary limit, Pst{V) = 
limt^ao P{V,t): while for smaller V the character of the original Levy stable behaviour is 
preserved (the original power-law behaviour, that is, persists to intermediately large V), 
for even larger V the corrections due to the dissipative non-linearity are visible in the 
transition(s) to steeper slope(s). 

These dissipative non-linearities remove the divergence of the kinetic energy from the 
measurable subsystem of the random walker. In the ideal mathematical language, the sur- 
rounding bath provides an infinite amount of energy through the Levy noise, and the cou- 
pling via the non-linear friction dissipates an infinite amount of energy into the bath, and 
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FIG. 10 Variance {V'^{t)) as function of time t, with the quartic term set to zero, 74 = and 
7o = 1.0 for all cases. The variance is finite for the cases /i = 2.0,72 = 0.0; ^ = 1.2,72 = 0.1; 
/i = 1.5,72 = 0.1; and fj, = 1.8,72 = 0.1. These correspond to the left ordinate. For the case 
fi = 1.2, 72 = 0.0, the variance diverges, strong fluctuations are visible; note the large values of this 
curve corresponding to the right ordinate. 
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thereby introduces a natural cutoff in the Icinetic energy distribution of the random walker 
subsystem. Physically, such divergencies are not expected, but correspond to the limiting 
procedure of large numbers in probability theory. We showed that both statements can be 
reconciled, and that Levy processes are indeed physical. 

Also Gaussian continuum diffusion exhibits non-physical features, possibly the most 
prominent being the infinite propagation speed inherent of the parabolic nature of the 
diffusion equation: even at very short times after system preparation in, say, a state 
P{x,0) = S{x), there has already arrived a finite portion of probability at large x. This 
problem can be corrected by changing from the diffusion to the Cattaneo (telegrapher's) 
equation. Still, for most purposes, the uncorrected diffusion equation is used. Similarly, 
one often uses natural boundary conditions even though the system under consideration is 
finite, since one might not be interested in the behaviour at times when a significant portion 
of probability has reached the boundaries. In a similar sense, we showed that "somewhere 
out in the wings" Levy flights are naturally cut off by dissipative non-linear effects. How- 
ever, instead of introducing artificial cutoffs, knowing that for all purposes Levy flights are 
a good quantitative description and therefore meaningful, we use "pure" Levy stable laws 
in physical models. 



H. Bi-fractional transport equations 



The coexistence of long-tailed forms for both jump length and waiting time PDFs was 
investigated within the CTRW approach in Ref. (11221 ). discussing in detail the laminar- 
localised phases in chaotic dynamics. In the framework of fractional transport equations, 
the combination of the waiting time PDF ip{t) ~ (0 < a < 1) and jump length PDF 

A(x) ~ leads to a d ynam ical equation with fractional derivatives in respect to both 

time and space (0; [sol; [giiTTlsh: 



d_ 
dt 



P{x,t) = Ki^\Dl- 



d\x\^ 



■P{x,t). 



(62) 



As long as the condition l<a</i<2is met, both exponents can be chosen within the 
entire range (jsol ). For yU = 2, in particular, this equation covers both sub- and superdiffusion 
up to ballistic motion, the latter corresponding to the wave equation (jssi ). A closed form 
solution of Eq. (1621) can be found in terms of Fox's if-functions, see Ref. (jo^ ). where also 
some special cases permitting elementary solutions are con s idere d. Bi-fractional diffusion 

T[iQi: 



equations were also discussed in Refs. (I32l : l35l : 160 
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Planck equation with a pow er-law dependence oc 
was studied in Refs. ( l55l : |74| ). 



\x\ 



). A bi-fractional Fokker- 
of the diffusion coefficient 



I. Levy walks 



Levy walks correspond to the spatiotemporally coupled version of continuous time random 
walks. The waiting time and jump length PDFs, that is, are no l ong er decoupled but appear 
as conditional in the form tp{x,t) = X{x)p{t\x) (or ijj{t)p{x\t)) (1671 ). In particular, through 
the coupling p{t\x) = ^6 {\x\ — vt"), one introduces a generalised velocity f , which penalises 
long jumps such that the overall process, the Levy wal k, att ains a finite variance and a PDF 
with two spiky fronts successively exploring space (69; 120l ). Thus, Levy walks have similar 
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properties to generalised Cattaneo/telegraphers' equation-type models (1491 : l82l : l89l ). As we 



here focus on the properties of transport processes governed by the Langevin equation under 
Levy noise, we only briefly introduce Levy walks. 

On the basis of fractional equations, formulations were obtained for the description of 
Levy walks in the presence of non-trivial external force fields, with the same restriction 



to lower order moments in respect to an Le vy w alk process (jSSi : 1871 ) . Recently, however. 



a coupled fractional equations was reported (11071 ). which describes a force-free Levy walk 
exactly. Thus, it was shown that the fractional version of the material derivative d/ dt±d/dx, 



,DfP(x±t,t), 



defined in Fourier-Laplace space through 



^ 1^ [dif{x, t);uy,k^ = {u± ikffik, u) 



(63) 



(64) 



acts on X and ^ on t) replaces the uncoupl ed fr actional time operators, see also the 
detaile d dis cussion of Levy walk processes in Ref. (Hl^). Although one may argue for certain 
forms f IO7I ). there is so far no derivation for the incorporation of general external force fields 
in the coupled formalism. We note that a very similar fractional approach to Levy walks 
was suggested in Ref. (18 ll ). 



IV. SUBDIFFUSION AND THE FRACTIONAL FOKKER-PLANCK EQUATION 

A. Physical foundation of subdiffusion 

The stochastic motion of a Brownian particle of mass m is described by the Langevin 
equation (0; 27; 73) 

dj^x dx 
m-— = —mrjiv + F{x) + mr{t), v = — , (65) 

where F{x) is an external force field, and 771 is the friction coefficient. The erroneous 
bombardment through the surrounding bath molecules is described by the fluctuating noise 
r{t). To properly describe Brownian motion, r[t) has to be chosen (5-correlated (white) 
and Gaussian distributed. This is, the time averages of r{t) are: r{t) = 0, and r(t)r{t') = 
DS{t — t'), where D is the noise strength. After averaging of the fluctuations, the velocity 
moments become ml 



{Av) =(t^v- ^\ At, {{Av)') = ^^^At + O ([Af ) . (66) 
\ m J m 

Both are proportional to At. For regular Brownian motion, these increments are used as 
expansion coefficients in the Chapman-Kolmogorov equation (jol). 

Subdiffusion now comes about through a so-called trapping scenario. Trapping describes 
the occasional immobilisation of the test particle for a waiting time distributed according 
to the distribution ip{t). Here we assume that the particle leaves the trap with the same 
velocity it had prior to immobilisation (this condition can be relaxed). In between trapping 
events, the particle is assumed to follow the regular Langevin equation such that each 
motion event on average lasts for a mean time r*. Choosing ?/'(t) with a finite characteristic 
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waiting time one can show that this trapping scenario indeed preserves Brownian motion. 
However, once the characteristic waiting time diverges, it can be shown by apphcation of the 
generalised central limit theorem that the occurrence of a large number of trapping events 
leads to subdiffusion and the fractional Klein-Kramers equation for the joint PDF P{x, v, t) 



(1861 : 192I : I93I ). Integrating out the velocity coordinate, the fractional Fokker-Planck equation 



emerges for the PDF P{x,t): 

In mathematical terms the trapping scenario corresponds to a subordination principle: 
Trapping events cause a transformation of the internal step time of the random walk to a 
different observation time, corresponding to the relation 

POO 

P{x,t)= S'ais,t)Pi{x,s)ds, (68) 
Jo 

where Pi{x,t) is the solution of the regular Brownian Fokker-Planck equation with a = 1. 
Explicit forms of the kernel (S'a{s,t) are known in terms of Fox' if-function (jlTI ). In the 
Laplace domain, the kernel (S'a{s,t) has the comparatively simple form 

^„(.,«)^^expf-^«".y (69) 

Note that this transformation guarantees the existence and positivity of the solution P{x, t) 
of the fractional Fokker-Planck equation if only the corresponding Brownian problem pos- 
sesses a well-defined solution. We note that recently an alternative approach t o co ntinuous 
time random walk subdiffusion in terms of (5-noise spike trains was presented f lOll ). 



B. Linear response and fluctuation-dissipation relation 



From the fractional Fokker-Planck equation (!67j) it is straightforward to prove that in the 
presence of a constant field Fq the linear response relation 



{x{t)) 



Fq 



-{x\t)) 



F=0 



(70) 



is fulfilled ( 11 71 : l83l ). This is due to the fact that the waiting time distribution %l)[t) is 
independent of the force. Given the comparatively large traps required to create the long- 
tailed form of ip{t), this assumption appears reasonable for not too large Fq. Experimentally, 
the linear response relation ( !70l) was verified (1611 ). 

The stationary solution of the fractional Fokker-Planck equation fl671) is the standard 
Boltzmann-Gibbs equilibrium, limt^oo -P(a^, ^) = =^exp{— 1^(x)/[A;b7']}, where ^ is a nor- 
malisation constant (1831 ). This result is immediately obvious from the picture of subordi- 
nation, that changes the temporal spacing of events but preserves the causal mode relax- 
ation. As F{x) = —V'{x), from the stationary solution of Eq. (I67j) by comparison with the 
Boltzmann-Gibbs form one recovers the generalised Einstein-Stokes relation (1831 ) 



mrjo 



(71) 



reflecting the preservation of the fluctuation-dissipation relation for the subdiffusion process. 
An experimental verification of this relation was reported by dlT 



25 




2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 






-1.5 -1 



-0.5 0.5 

X 



1 1.5 



FIG. 12 Propagator P{x,t) of the fractional diffusion equation with a = 1/2 (left) and the normal 
diffusion equation (right) for consecutive times. 



C. Propagator 

In absence of an external force, the fractional diffusion equation is solved by Fox' H- 
f unction (1171 ) . 
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(72) 



an equivalent form to the original solution by (11051 ). Its series expansion reads (1171 ) 

„2 \ ri/2 
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and asymptotically reaches the stretched Gaussian form 

1 



(73) 



W{x,t) 



1 



2- a /a\"/(2-") 
Y~ V2 



2\ (l-a)/(2-a) 



-(l-a)/(2-a) 



X exp 



\x\ 


l/(l~a/2)\ 







(74) 



valid for |x| ^ y/Kat"'. The l atter corresponds to the known result from continuous time 
random walk theory (1691 : Il20l ). The if- function simplifies if the exponent a is a rational 
number. For instance, for a = 1/2, it can be rewritten in terms of the Meijer G-function 
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(75) 



that is known in symbolic mathematics packages such as Mathematica. Using this repre- 
sentation in Fig. [12] we show the subdiffusive propagator with its pronounced cusps at the 
site of the initial condition, in comparison to the smooth Gaussian propagator of Brownian 
motion. 
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D. Boundary value and first passage time problems 

In the presence of reflecting, absorbing, or mixed, boundary conditions, the narrow jump 
length distribution of subdiffusive processes makes it possible to use analogous techniques to 
calculate the propagator as known from normal diffusion. Thus, separation into eigenmodes 
or the method of images can be applied, the only difference entering through the time- 
dependence of the modes. A very convenient way is to employ the Brownian result for 
a given geometry, and subordinate that result to obtain the behaviour for a subdiffusing 
particle, according to Eq. (168|) . 

In the presence of an absorbing boundary an important quantity characterising the dy- 
namics of the system is the survival probability 

y^{t) = I ^{x,t)dx, (76) 



where D denotes the interval over which ^(x, t) is defined. The initial value of the survival 
probability is ^^(0) = 1, and it decays to zero for long times. From ^a{t), we can define 
the first passage time density 

Pait) = -J^^it). (77) 

The mean first passage time is 

I'OO 

^ = Pa{t)tdt. (78) 

Note that from Eqs. (177|) and (l68l) . we obtain the relation 

Pa(n)=pi(n") (79) 

between the subdiffusive and Brownian results in Laplace space. Thus, these are connected 
by a relation similar to the subordination (1681) . with kernel exp (— sm") instead of S'a{s,u) 



(ll8|;l88D. 

For the three most prominent cases of first passage time problems, we obtain the following 
subdiffusive generalisations: 

(i) For subdiffusion in the semi-infinite domain with an absorbing wall at the origin and 
initial condition P{x, 0) = 5{x — xq) it was found that (jool ) 

\V[-a/2)\Kc^ 

i.e., the decay becomes a flatter power-law than in the Markovian case where p{t) ~ 

(ii) Subdiffusion in the semi-infinite domain in th e presence of an external bias V falls 



off faster, but still in power-law manner (1341 : l88l : |10J): 



p{t) ~ (81) 

In strong contrast to the biased Brownian case, we now end up with a process whose char- 
acteristic time scale diverges. This is exactly the mirror of the multiple trapping model, 
i.e., the classical motion events become repeatedly interrupted such that the immobilisation 
time dominates the process. In contrast, for a = 1 the result 

m / (xq -Vtf \ . 
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is valid, producing the classical form T = Xq/V for the mean first passage time, 
(iii) Subdiffusion in a finite box (jool ): 



p{t) ~ t 



-l-o 



(83) 



i.e., this process leads to the same scaling behaviour for longer times as found for the biased 
semi-infinite case (ii). 

The latter two results should be compared to the classical Scher-Montroll finding for 
the first passage time density of biased motion in a finite system of size L with absorbing 
boundary condition. In that case, the first passage time density exhibits two power-laws 



p{t) 



r-\ t < T 

t-^-", t > T 



(84) 



the sum of whose exponents equals —2 (99; 9^; 103 : 104 ). Here, r is a system size dependent 
time scale (1991). 



E. Fractional Ornstein-Uhlenbeck process 

The Ornstein-Uhlenbeck process corresponds to the motion in an harmonic potential 
V{x) = |mu;^x^ giving rise to the restoring force field F{x) = —mu^x, i.e., to the dynamical 
equation 

|p(.,*)^„A-(|;^ + if.^)p(.,*). (85) 

From separation of variables, and the definition of the Hermite polynomials ([l]), one finds 
the series solution for the fractional Fokker-Planck equation with the Ornstein-Uhlenbeck 



potential (1171 : 1831 ) 



mux \ I muj'^x'^ 



plotted in Fig. [131 Individual spatial eigenmodes follow the ordinary Hermite polynomials of 
increasing order, while their temporal relaxation is of Mittag-Leffler form, with decreasing 
internal time scale {'r]a/[nuj'^])^^" . Numerically, the solution (l86l) is somewhat cumbersome 
to treat. In order to plot the PDF P{x,t) in Fig. [131 it is preferable to use the closed form 
solution (we use dimensionless variables) 

P(x, t) = .^^^ exp I - - , — 1 (87) 



^27r(l - 




of the Brownian case, and the transformation ( 1681) to construct the fractional analogue. 

Fig. [131 shows the distinct cusps at the position of the initial condition at Xq = 1. The re- 
laxation to the final Gaussian Gibbs-Boltzmann PDF can be seen from the sequence of three 
consecutive times. Only at stationarity, the cusp gives way to the smooth Gaussian shape of 
the equilibrium PDF. By adding an additional linear drift V to the harmonic restoring force. 
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FIG. 13 Time evolution of the PDF of the fractional Ornstein-Uhlenbeck process {a = 1/2). The 
initial condition was chosen as 6{x — 1). Note the strongly persistent cusp at the location of the 
initial peak. Dimensionless times: 0.02, 0.2, 2. The dashed line corresponds to the Boltzmann 
equilibrium. 

the drift term in the FFPE (1671) changes to —d{x — V)P{x, t)/dx, and the exponential in ex- 
pression (1871) takes on the form exp {—[x — V — (xq — V)e~^] / [2 (1 — e"^*)]). As displayed 
in Fig. HH the strong persistence of the initial condition causes a highly asymmetric shape 
of the PDF, whereas the Brownian solution shown in dashed lines retains its symmetric 
Gaussian profile. 

Let us finally address the moments of the fractional Ornstein-Uhlenbeck process, Eq. fl86l) . 
These can be readily obtained either from the Brownian result with the integral transforma- 
tion ( l68l) . or from integration J dxx^- of the FFPE ( l67l) . For the first and second moments 
one obtains: 



respectively. The first moment starts off at the initial position, xq, and then falls off in a 
Mittag-Leffler pattern, reaching the terminal inverse power-law ~ t~". The second moment 
turns from the initial value Xq to the thermal value x"^^ = ksT / {muj'^). In the special case 
Xq = 0, the second moment measures initial force- free diffusion due to the initial exploration 
of the fiat apex of the potential. We graph the two moments in Fig. [T5] in comparison to 
their Brownian counterparts. 

V. FUTURE DIRECTIONS 




(88) 



and 




(89) 



Anomalous diffusion is becoming widely recognised in a variety of fields. Apart from the 
anomalous spreading of tracers and the consequences for the propagator, additional questions 
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P (x,t) 




FIG. 14 Time evolution of the PDF of the fractional Ornstein-Uhlenbeck process with superposed 
constant force of dimensionless strength V = —4 (a = 1/2). The initial condition was chosen as 
(5(x — 4). Dimensionless times: 0.02, 0.2, 2. The dashed lines corresponds to the Brownian solution 
at times 0.5 and 50 (in essence, the stationary state). Again, note the cusps due to the initial 
condition, causing a strongly asymmetric shape of the PDF in contrast to the Gaussian nature of 
the Brownian counterpart. 



such as ageing and weak ergodicity breaking become important for the understanding of 
experiments and their modelling on complex systems. At the same time, the theory of 
anomalous processes is expanding. For instance, regarding questions on the weak ergodicity 
breaking, the calculation of multipoint moments or the correct introduction of cutoffs are 
currently being worked on, to complete the world Beyond Brownian Motion (1681 ). 
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FIG. 15 First (xq = 2, dashed line) and second {xq = 0, full line) moment of the fractional 
Ornstein-Uhlenbeck process (a = 1/2), in comparison to the Brownian case. log^Q— log^^Q scale. 
The dotted straight lines show the initial (sub)diffusive behaviour with slopes 1/2 and 1, in the 
special case xq = chosen for the second moment. 
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